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Abstract - We study mesoscopic fluctuations in a system in which there is a continuous connection 
between two distinct Fermi liquids, asking whether the mesoscopic variation in the two limits is 
correlated. The particular system studied is an Anderson impurity coupled to a finite mesoscopic 
reservoir described by random matrix theory, a structure which can be realized using quantum 
dots. We use the slave boson mean field approach to connect the levels of the uncoupled system 
to those of the strong coupling Nozieres Fermi liquid. We find strong but not complete correlation 
between the mesoscopic properties in the two limits and several universal features. 



Introduction. — The Fermi liquid is a ubiquitous 
state of electronic matter [lr[3]. Indeed, it is so com- 
mon that systems can have several different Fermi liquid 
phases in different parameter regimes (controlled by dif- 
ferent fixed points), leading to cross-overs between Fermi 
liquids with different characteristics. Examples of such 
cross-overs include, for instance, the half-filled Landau 
level (high-temperature to low-temperature connection) 
[J, heavy fermion materials OIMT], as well as the sim- 
ple spin 1/2 Kondo problem which will be our main con- 
cern in this paper. In the bulk, clean case, the evolution 
of the quasi-particles in such a cross-over is straight for- 
ward: both sets of quasi-particles are labeled by k be- 
cause of translational invariance and so are in one-to- 
one correspondence. However, in the absence of trans- 
lational invariance — such as in a disordered or mesoscopic 
setting — interference affects the two sets of quasi-particles 
differently. In such a situation, it is interesting to ask how 
the quasi-particles in one Fermi liquid are related to those 
in the other. 

The Kondo problem provides a particularly clear ex- 
ample: at weak coupling (high temperature) the elec- 
trons in the Fermi sea are nearly non-interacting while the 
strong coupling (low temperature) behavior is described 
by Noziere's Fermi liquid theory [HE]. The connection 
between high and low temperature is provided, e.g., by 
Wilson's renormalization group calculation [9], yielding a 
smooth cross-over. 



We now break translational invariance by supposing 
that the size of the Fermi sea is finite; it could consist 
of, for instance, a large quantum dot or metallic nano- 
particlc. The density of states in the electron sea will 
typically have low energy structure and features, in con- 
trast to the intensively studied flat band case. The finite 
size effects introduce two additional energy scales: (i) a 
finite mean level spacing, leading to what is called the 
"Kondo box" problem [TMT2"] . and (ii) the Thouless en- 
ergy £/rh = h/r c where r c is the typical time to travel 
across the finite reservoir. When probed with an energy 
resolution smaller than £/rh, both the spectrum and the 
wave- functions of the electron sea display mesoscopic fluc- 
tuations, which affect the Kondo physics p~3VH5J . Disorder 
in the electron sea causes similar effects fTBHIIj] . 

Consider the system shown in Fig. [1] a small Kondo 
dot coupled to a large "reservoir dot" probed weakly by 
tunneling from a tip. In the high temperature regime, 
the small dot is weakly coupled to the large dot which 
is essentially non-interacting. Mesoscopic fluctuations of 
the density of states translate into mesoscopic fluctua- 
tions of the Kondo temperature. Once this translation 
is taken into account, the high-temperature physics re- 
mains essentially the same as in the flat band case [T3HT5] ; 
in particular, physical properties can be written as the 
same universal function of the ratio T/Tk as in the bulk 
flat-band CMSC . MS long as Ik is understood as a realiza- 
tion dependent parameter |13H15| . In contrast, the con- 



p-1 



Dong E. Liu et al. 



Weak Coupling : T»T k 

e 



Strong Coupling : T«T k 




Non interacting electron 



Non interacting electron 
at SC fixed point 



Fig. 1: (Color online) Schematic illustration of small-large 
quantum dot system. Left panel : weak coupling limit T 3> Tk. 
Right panel : strong coupling limit T <C Tk. The energy levels 
and wave functions probed by the tip change from one Fermi 
liquid regime to the other. 



sequences of mesoscopic fluctuations for low temperature 
Kondo physics (T <C Tk, strong coupling limit) remain 
largely unexplored. A few things are nevertheless known: 
for instance, the very low temperature regime should be 
described by a Nozieres-Landau Fermi liquid, as in the 
original Kondo problem. Indeed, the physical reasoning 
behind the emergence of Fermi liquid behavior at low tem- 
perature, namely that for energies much lower than Tk 
the impurity spin has to be completely screened, applies 
as well in the mesoscopic case as long as T < A <c Tk 
[20 2S]. In this case, the system consists of the Kondo 
singlet plus non-interacting electrons with a 7r/2 phase 
shift as shown in the right panel of Fig. [T] 

Measurements of the conductance through the large dot 
or the ac response to the tip reveal the mesoscopic fluc- 
tuations of the energy levels and wavcf unctions [26 | l27 j . 
Thus, such experiments can probe the connection between 
the quasi-particles in the two Fermi liquid regimes, as well 
as the properties of the intermediate strongly correlated 
Kondo cloud [2"2Tl2"5] . In this paper, we study this con- 
nection explicitly, using slave boson mean field (SBMF) 
theory [51IMTl l2"51 - f52"] to treat the interactions and random 
matrix theory (RMT) [321 to model the mesoscopic fluc- 
tuations. We find that the correlation between the prop- 
erties of the two sets of quasi-particles is substantial but 
not complete. 



Model. — The system pictured in Fig. [T] can be de- 
scribed by the Hamiltonian H = -Hbath + iTmp where 
-ffbath describes the mesoscopic electronic bath and H lmv 
describes the local "magnetic impurity" — small quantum 
dot, nanoparticle, or magnetic ion — and its interaction 
with the bath. i?bath is the non-interacting Hamiltonian 
#bath = J2i,a( e i - mKoAo-' where i = 1, • • • ,N labels the 
levels, a =t, l is the spin component, and \i is the chemical 
potential. H lmp is 

H imp = V E[ c o<A + 4coa] + E d E dido (1) 



where the d a operators refer to the impurity site with en- 
ergy Ed and the position of the impurity is taken to be 
r = 0. We take the local Coulomb interaction between d- 
electrons to be U — oo; thus, states with two c?-electrons 
must be projected out. Finally, the local electronic op- 
erator co<j is related to the bath eigenstate operators c ia 
through c 0(T = Y,i=i 01(0)^ where 4 (r) = (r|i) are the 
one-body wave functions of i?bath with the local normal- 
ization relation J^i I §i (0)1 2 = 1- 

To study the mesoscopic fluctuations, we assume that 
the classical dynamics within the large dot is chaotic, and 
thus that the energy levels and the wave functions at 
the impurity site <fo(0) are described by random matrix 
theory (RMT) [15l[Ml[33j, specifically, by the Gaussian 
orthogonal ensemble (GOE) for time reversal symmetric 
systems and the Gaussian unitary ensemble (GUE) for 
systems in which time reversal is broken [331I34]. 

Applying the SBMF approximation jMH71l2"5H52] . we in- 
troduce auxiliary boson b* and fermion j\ operators, such 
that d a = b'fer, with the constraint w b + /J/^ = 1. 
Since the Hamiltonian is invariant with respect to a U(l) 
gauge transformation, the bosonic field can be treated as 
a real number : b, \— > rj. The constraint condition is sat- 
isfied by introducing a static Lagrange multiplier, £. One 
thus obtains the SBMF effective Hamiltonian 
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The mean field parameters r/ and £ are obtained by min- 
imizing the free energy of the system, taking fi = 0. Us- 
ing the equations of motion from the mean-field Hamilto- 
nian Eq. ([2]), we obtain, after some algebra, the imaginary 
time Green function 
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from which all the properties of the system can be derived. 

The eigenvalues A K and eigenstates \ip K ) (k = 
0, 1, • • • , N) of the mean field Hamiltonian Eq. (|2|) corre- 
spond to the quasi-particles of the strong coupling limit. 
Because the low temperature/energy regime of the sys- 
tem is a Fermi liquid, the mean field approach provides a 
good description of the low energy properties of the strong 
coupling limit, but it is not expected to be accurate at 
higher energies. As a consequence, it is mainly the range 
|A K — fi\ < Tk which is physically relevant in terms of 
Kondo physics. We shall therefore in the following con- 
centrate on this energy range. Since the tunneling strength 
at energy E between an external tip and the large quan- 
tum dot (see Fig. [TJ depends on the line-up of the levels 
and the wavefunction intensity, both A K and |'0K(r)| 2 are 
measurable in experiments. 
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We now study the relation between the {X K ,\ip K )} and 
the {ei,\4>i)}. Expressing the Green function of Hmf as 
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one sees that (A K — fi) are the poles of the Green function 
Gff(z) = (f\G(z)\f) . Eq. ((3j) then immediately implies 
that the A K are solutions of the equations 
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where £o(£) = Ed+p — £ (interpreted as the position of the 
Kondo resonance if the system is in the Kondo regime) and 
r c ff = ttpo7] 2 Vq (interpreted as the width of the Kondo 
resonance, which gives the scale of the Kondo tempera- 
ture), po = 1/A is the mean density of states. Note that 
Eq. ([S]) implies that there is one and only one X K in each 
interval [ei,e i+1 ]. 

The probability of overlap between the eigenstate k and 
the impurity state |/}, u K = |(/|^ K )| 2 , is a key ingredient 
in how the wave function amplitude at r is affected by the 
Kondo singlet. Since the u K are the residues of Gg(z), 
Eq. © implies 
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For |A — £o(£) ^ T e ff, one contribution dominates the 
sum on the left hand side of Eq. (J5j — namely, the closest 
€i to A K , call it i(n) — in which case X K = £j( K j + <5 K A with 
S K ~ r cff |(/) 4(K) (0)| 2 /[7r(A K - So)} < 1- As expected, the 
two spectra nearly coincide. Similarly, the participation 
of the wavefunctions in the singlet state is small: from 
Eq. © 

r off l<W0)| 2 A A 
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In contrast, for | A — (0 1 ^ FeS, the right hand side of 
Eq. ([5J can be neglected. The typical distance between a 
A K and the closest ti is then of order A, and u K ~ A/r c ff. 
In the limit Ik ~ r c ff ^ A, only the wave function 
amplitudes for energy levels within the Kondo resonance, 
|A K — £o(£)l ^ r o ff, will be significantly affected. 

Energy Spectral Correlation. — To characterize 
the relation between the weak and strong coupling levels, 
{e^} and {A K } respectively, we consider the distribution of 
the normalized level shift defined by 
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where and e^+i are the two levels which sandwich A K . 
The range of 5" is from to 1. The probability distribution 
P(S) obtained numerically using the SBMF approxima- 
tion by sampling a large number of realizations is shown 
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Fig. 2: (Color online) The distribution of S, including both 
\X K - ti\/\e l+1 - a\ and |A K - e i+ i|/|e i+ i - q\, from the SBMF 
treatment of the infinite-(7 Anderson model; (a) GOE, (b) 
GUE. The dashed lines are the results of the toy model. In- 
set (a): the cumulative distribution F(S) = J Q S p(x)dx of the 
Vo — 1-0 GOE data compared to the toy model. Note the 
presence of the square-root singularity. Parameters: full band 
width D — 3, impurity energy level Ed = —0.7, 500 energy 
levels within the band, and 5000 realizations used. 



in Fig. [3] for several cases. Only the levels that are within 
the Kondo resonance are included; that is, levels satisfying 
\X K — £q\ < r o ff/2. Note in particular two features of the 
numerical results: (i) the strong coupling levels are more 
concentrated near the original levels in the case of the 
GOE while they are pushed away from the original levels 
in the GUE, and (ii) the distribution found is completely 
independent of Vq . 

An explanation for both of these features can be found 
from a simple analytic approximation to the distribution 
P(S). Well within the resonance, \X K — £q\ <S r c ff, the 
r.h.s. of Eq. ([5]) can be set equal to zero, thus leading to 
the simplification J2i=i l ( M0)| 2 /(A K — (-i) ~ 0. Focusing 
on the level A K located between and e;+i, we consider 
a toy model in which the influence of all but these closest 
e's is neglected, yielding the much simpler equation for A K 
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In RMT, the wavefunction amplitudes \(j>i\ 2 and |<^i+i| 2 
are uncorrclated and distributed according to the Porter- 
Thomas distribution I331I34| . Notice that all energy scales 
(Vq, A, etc.) have disappeared from the problem except 
for Se = e.; + i — ej. The resulting distribution of A K is there- 
fore universal, depending only on the symmetry under 
time reversal. Hence the empirical observation in Fig. [5] 
that the curves are independent of Vq . 

Integration over the Porter-Thomas distributions gives 
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Breaking time-reversal symmetry thus affects drastically 
the correlation between the low temperature levels A K and 
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the neighboring high temperature ones ej and £j+i. Time- 
reversal symmetric systems show clustering, with a square 
root singularity, of the A K 's close to the ej's, while for sys- 
tems without time-reversal symmetry the distribution is 
uniform between ej and The GUE result can be im- 
proved by taking into account the other levels on average; 
this yields the expression plotted in Fig. [U[b), but as it is 
lengthy we do not specify it here. Note that this improved 
toy model does give the bunching of levels in the middle 
of the interval seen in the numerics. 

The difference between P(S) in the two ensembles 
comes from the very different wave function distribution: 
the GOE Porter- Thomas distribution has a square root 
singularity at 0i(O)| 2 = while it is finite for the GUE. 
The high probability of small wave function amplitudes in 
the GOE leads to the clustering of strong coupling levels 
around the original ones. To explore this in the SBMF nu- 
merical results, we plot the cumulative distribution func- 
tion on a log-log scale in the inset in Fig. [2] the resulting 
straight line parallel to the toy model result shows that, 
indeed, the square root singularity is present. 

Wave Function Correlations. — A key quantity in 
quantum dot physics is the magnitude of the wave func- 
tion of a level at a point in the dot that is coupled to an 
external lead (see Fig. [T]). This quantity is directly re- 
lated to the conductance into the dot when the chemical 
potential in the lead is close to the energy of the level. 
We assume that the probing lead is very weakly coupled, 
so that the relevant quantity is the wave function in the 
absence of leads. To see how the tunneling to an outside 
lead at r is affected by the coupling to the impurity, we 
study the correlation between the strong coupling wave- 
function intensity |?/> K (i)(r)| 2 and its weak coupling coun- 
terpart |</)i(r)| 2 , with k(i) = i for A K < £ and = (i + 1) 
for X K > £q. Specifically, we consider the correlator 
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The average (•) here is over all realizations, for arbitrary 
fixed r / 0, and a(-) is the square root of the variance of 
the corresponding quantity. 

Results for C^ K u\ from the SBMF approach to the 
infinite-C/ Anderson model are shown in Fig. [3] Two ways 
of showing the dependence on the argument i are used: in 
the left hand panels, the x-axis is simply the (average) en- 
ergy from the middle of the band, namely Sli = iA — D/2, 
while in the right hand panels, this energy is scaled so that 
the x-axis is the energy from the center of the Kondo res- 
onance in units of the Kondo temperature (see caption for 
exact expression). Because the infinite-C/ Anderson model 
is inherently not particle-hole symmetric, the location of 
the Kondo resonance is not at zero but rather increases as 
V increases so that the average occupation of the impurity 
level is less than one. 

The scaled curves have a very natural interpretation. 
First, those states which do not participate in the Kondo 
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Fig. 3: (Color online) Wave function correlation, Ci jK u\, for the 
SBFM approach to the infinite-17 Anderson model, (a) GOE 
and (c) GUE, as a function of the average distance from the 
middle of the band, (b) GOE and (d) GUE, as a function of the 
rescaled average distance [iA — D/2 — (£o(£) — /z)]/r o ff. Dashed 
line: analytic approximation, Eq. (|16p . Parameters: full band 
width D — 3, impurity energy level Ed = —0.7, 500 energy 
levels within the band, and 2000 realizations. 



singlet state at low temperature, \5e — (Ed — £)| 3> r e g, 
are essentially unchanged, C i K ^ ~ 1. In contrast, those 
states with energies within the Kondo resonance are sub- 
stantially changed by interaction with the impurity. The 
universality of the low energy Kondo physics is nicely 
demonstrated by the collapse of all the numerical curves 
for different coupling strengths onto universal curves, one 
for the GOE and one for the GUE. 

The most interesting feature in Fig. [3] is that the corre- 
lation does not go to zero at the center of the Kondo res- 
onance, even for strong bare coupling. Clearly, the wave 
functions of the weak coupling and strong coupling Fermi 
liquid states are similar to each other in that the inter- 
ference pattern in the original wave function is not com- 
pletely wiped out by the formation of the Kondo singlet 
state. This residual correlation should be observable as a 
correlation in the conductance probed by an external tip. 

Expressing the wave function probability ^(r)! 2 as the 
residue of G(r, r), and using that in the semiclassical limit 
the magnitude of the unperturbed wave- functions <pi(r) 
are uncorrelated at different points and with the energy 
levels, one can show [35] that in the limit r o fj 3> A, the 
correlator is given by 



(13) 



where m = 7]Vo4>i(0). An approximation to Cj )K (j) can be 
obtained by taking the energy levels to be evenly spaced 
and replacing the wave function intensities by the average 
value; this yields 
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where S K = (X K (i) — £i)/A. Within the same approxima- 
tions, Eq. ([5]) then implies 

^% g ° = - ^ , 1 . = cotan(7r5 K ) . (15) 

Toff 7T ^ d R - J 

By defining A K = X K — £q, we obtain for the correlator 
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which, as anticipated, depends only on the ratio (A K /r e ff). 
The curve resulting from this expression is shown in Fig. [3] 
(b) and (d); it yields the value C itK ^ ~ 4/ it 2 at the mini- 
mum, independent of all parameters |35j . The value found 
numerically is slightly smaller but in reasonable agree- 
ment. 

Conclusion. — We have presented the first study of 
mcsoscopic fluctuations in two distinct but continuously 
connected Fermi liquids by using the slave boson mean 
field approximation to calculate the strong-coupling lev- 
els. In the specific case that we study — a small quantum 
dot coupled to a large reservoir quantum dot with chaotic 
dynamics — the fluctuations of single particle properties in 
the two limits are highly correlated, universal, and very 
sensitive to time reversal symmetry. Indeed, each strong 
coupling level must lie between two of the original lev- 
els (the spectra are interleaved), and, in the GOE case 
but not the GUE, the levels within the Kondo resonance 
are clustered about the weak coupling ones. Similarly, 
while the wave function correlation dips within the Kondo 
resonance, it remains substantial, showing that while the 
corresponding wave function is strongly affected by the 
Kondo screening it retains a surprisingly substantial over- 
lap with the original wave function. We expect a similar 
strong correlation between the properties of continuously 
connected distinct Fermi liquids in other systems. An in- 
teresting extension would be to study such correlations 
when a quantum phase transition intervenes as in, e.g., 
the two impurity Kondo model. 
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